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ABSTRACT 

We estimate the gravitational wave (GW) background from Population III (Pop III) stars using the results 
from our hydrodynamic simulations (Suwa et al. 2007). We calculate gravitational waveforms from matter 
motions and anisotropic neutrino emissions for single Pop III stars. We find that the GW amplitudes from 
matter motions are dominant until just after bounce, but those from neutrinos dominate later on at frequencies 
below ~ 10 Hz in the GW spectrum. Computing the overall signal produced by the ensemble of such Pop III 
stars, we find that the resultant density parameter of the GW background peaks at the amplitude of ficw ~ 10~ 10 
in the frequency interval ~ 1 - 10 Hz. We show that such signals, depending on the formation rate of Pop III 
stars, can be within the detection limits of future planned interferometers such as DECIGO and BBO in the 
frequency interval of ~ 0.1 — 1 Hz. Our results suggest that the detection of the GW background from Pop III 
stars can be an important tool to supply the information about the star formation history in the early universe. 
Subject headings: gravitational waves — supernovae: general — black hole physics — neutrinos — methods: 
numerical — hydrodynamics 
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1. INTRODUCTION 

The observation of gravitational waves (GWs) is one of the 
most important tasks for exploring the less well known parts 
of the Universe. Already several ground-based laser interfer- 
ometers (TAMA300, LIGO, and GEO600) are now operating 
and taking data in the frequency range of 10 Hz - 10 kHz, 
where rapidly-collapsing stellar objects accompanied with 
strong gravity such as neutron stars (NSs) and black holes 
(BH s), are one of the most promis ing sources of GWs (see, 
e.g.. lNewll200H iKotake et al.ll2006l for reviews). The Laser 
Interferometer Space Anntena (LISA), covering 10~ 4 -10~ 2 
Hz, will be launched in the near fu ture, and, more over, future 
space missions such as DECIGO dSeto et alj|2001l) and BBO 
( Unga relli et al.1120051) . targeting the ~ 0. 1 Hz regime, are be- 
ing planned. Possible targets of these experiments in the low 
frequency region are gravitational wave backgrounds (GWBs) 
from both astrophysical and cosmological origins. In particu- 
lar the low-frequency window is thought to be indispensable 
for the detection of the primordial G WB generated during the 
inflationary epoch (Maggiore 2000). 

In addition, various possible astrophysical sources of 
GWBs have been i nvestigated. Among them are core- 
collapse supernovae dde Arauio et a l. 2004), inspiral and co- 
alescence of compact binaries dFarmer & Phinneyll2003l). ro- 
tating NSs with non-axisy mmetric deformations dFerrari et all 
1999), gamma-ray bursts (Hiramatsu et al. 2005), and NSs 
with phase transitions dSigll 120061) (see references therein). 
Recently, iBuonanno et al.l d2005l) pointed out that the contri- 
bution from Population III (hereafter Pop III) stars is most im- 
portant because it could be as large as the inflationary GWB 
around . 1 Hz where DECIG O and BBO are most sensitive 
(see also lSandick et al.ll2006l) . 
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Here, Pop III stars are the first stars, which formed out of 
a pristine, meta l free gas produced by the big bang ( for re- 
views, see, e.g. iBromm & Larsonll2004t lGloverll2005h . Re- 
cent simulations indicate that the initial mass function (IMF) 
of Pop III stars is expected to have been dominated by very 
massive stars with mass lOOM^ (IN akamura & Umemura 
1 2001b lAbel etai]|2002t IBromm et alj|2002t lOmukai & Palla 
120031) . These very massive stars encounter the electron- 
positron pair creation instability after central carbon burn- 
ing, which reduces the thermal energy and induces gravi- 
tational collapse. For stellar mass of less than ~ 260M Q , 
rapid nuclear burning releases large energy sufficient to en- 
tirely disrupt the star as pair-instability supernova. More mas- 
sive stars, which also encounter pair-instability, are so tightly 
bound that the fusion of oxygen is unable to reverse inf all - 
Such stars are th ought to collapse to BHs dBond et al.lll984t 
iFryer et afll2001l) . which we pay attention to in this letter. So 
far, there have been only a few hydrodynamic simulations 
studying the gravitational collapse of the BH forming Pop 
III stars dFrver et alJl2001b iNakazato etaT] l2006: Su wa et all 
120071 iLiuet al.H2007h . 

Due to the l ack o f Pop III star collapse simulations, 
IBuonanno et al.l d2005l) had to calculate the GWB with an as- 
sumption that the GW spectrum of a single Pop III star has 
the same shape as the one from an ordi nary core-collapse su- 
pernova. They employed the result of Miill er et~aT] d2004l) as 
the template. In order to determine the normalization of the 
GW amplitudes, they simply amplified the GW amplitudes of 
th e supernova by a f actor of ~ 1000, relying on the results 
by lFryer et alJ d2001). Needless to say, these assumptions and 
treatments can be validated by a consistent treatment, namely 
by hydrodynamic collapse simulations of Pop III star with the 
analysis of the GW wavefo rms. As for the mass an d forma- 
tion history of Pop III stars, IBuonanno et al.l d2005l) assumed 
that all stars formed with the same mass of 300M Q and at 
th e same redshif t as z ~ 15. This assumption is improved 
by lSandickeTaTI d2006l) . who calculated GWB from Pop III 
stars with the cosmic star formation history in t he framework 
of hier archical structure formation. However, ISandick et all 
(2006) also employed t he sam e spectrum of a single Pop III 
star as Buo nanno et all (12005b . 
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The purpose of this letter is to give an estimate of the GWB 
spectrum from Pop III stars by calculating the GW waveforms 
based on the results of hydrodynamical core-collapse simula- 
tions of Pop III stars. This letter is organized as follows: In 
the next section the numerical model is briefly described. In 
§3, we calculate the gravitational wave signal and its spec- 
trum of a single Pop III star collapse. In §4, we present the 
numerical result of the gravitational wave background from 
Pop III stars. §5 is devoted to summary and discussion. 

2. METHOD 

The numerical met hods are basically the same as the ones 
in our previous paper (Suwa et al. 2007). With the ZEUS-2D 
code (IStone & Normanll 19921) as a base for the hydro solver, 
we employ a realisti c equation of s tate b ased on relativistic 
mean field theory (IShen et al .11 19981) (see iKotake et al.ll2003l 
for the implementation) and treat the neutrino cooling by a 
multiflavour leakage scheme, in which 6 spe cies of neutrinos 
with pair, photo, and plasma processes by lltoh et all (1 19891) 
in addition to the standard charged current neutrino cooling 
reactions are included. Spherical coordinates (r, 8) with log- 
arithmic zoning in the radial direction and regular zoning in 
8 are used. One quadrant of the meridian section is covered 
with 300 (r)x 30 (8) mesh points. In our 2D calculations, 
axial symmetry and reflection symmetry across the equato- 
rial plane are assumed. We also calculated some models with 
60 angular mesh points, but found no significant differences. 
Therefore, we will report in the following the results obtained 
from the models with 30 angular mesh points. 

The initi a l co ndition is provided in the same manner of 
ISuwa et ail (120071) . We produced hydrostatic cores of 300, 
500, 700, and 900 M@, with the assumpti on of isentropic , 
whose values of entropy are taken from iBond et alj d!984l) . 
and with the constant electron fraction of Y e = 0.5. The sup- 
posed rotation law is cylindrical rotation, with the strength of 
rotational energy taken to be 0.5% of the gravitational energy 
in all models. 

In this study, we estimate the gravitational wave emis- 
sion from aspherica l mass motions via the Newtonian 
quadrupole formula of Monchmeve r et all (1 199 lb . In addition, 
we compute the gravitational wave strain from anisotropic 
neutrino em ission employing th e formalism introduce d by 
Epstein! (1 19781) and developed by Miill er & Jankal d 1997b and 
Kota ke et all d2007l) . The GW emission from neutrinos is 
given as 

Dh TT (t)=— a(t')L u (t')dt', (1) 

c J -oo 

where D is the distance to the source, h TT is the transverse- 
traceless and dimensionless metric strain, G is the gravita- 
tional constant, c is the speed of light, a(t) is the time- 
dependent neutrino-anisotropy parameter, and L v (t) is the to- 
tal neutrino luminosity. For the estimation of a, neutrinos are 
assumed to be emitted radially in each angular bin. 

3. GRAVITATIONAL WAVE OF A SINGLE POP III STAR COLLAPSE 

In this section, we discuss the GW emission from a single 
3OOM0 Pop III star collapse. FigureQ]depicts the strain versus 
time after bounce for a 300Mq star. The matter contribution 
dominates h TT during first 10 ms after bounce. Afterward, the 
neutrino part begins to contribute because the thermal shock 
occurs and the neutrino luminosity increases in the hot region 
behind the shock wave. At about 70 msec after bounce, the 
neutrino contribution converges to a constant value and the 
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FIG. 1 . — The gravitational wave strain, h TT , times the distance to the Pop 
III star, D, versus time after bounce (in seconds). Shown are the neutrino 
(blue dotted) and matter motion (green dashed) components along with their 
sum (red solid). 

matter contribution goes to zero. This represents the epoch of 
the BH formation. After BH formation, we do not calculate 
the GW emission from matter but dump with the timescale of 
light crossing time ~ 0(Rbh/c), where Rbh is the radius of the 
BH. This procedure does not affect the discussion of the fol- 
lowing section because the matter contribution of total GW is 
only in the high frequency region, which is not the main point 
of this letter. The total energy emitted in gravitatio nal waves is 
~2 x 10 50 erg, which is smaller than the result of lFrver et alj 
d2001l) by a factor of 10. The discrepancy seems to come from 
the difference of initial angular momentum distribution. 

Now we move on to discuss the features of the waveforms 
by the spectral analysis. From Figure|2] one can see the domi- 
nance of the neutrino-oriented gravitational waves at frequen- 
cies below ~ 10 Hz. This is because the waveform from neu- 
trinos shows the slower time-variability in comparison with 
the one from matter motion (see FigureQ]). On the other hand, 
the contribution of matter motion dominates the spectrum in 
the higher frequency region. The peak amplitudes from mat- 
ter motion and the frequencies of 10s of Hz are roughly at 
the same characteristic amplitude and frequency produced by 
the collapse o f the star to form a BH, which are evaluated in 
I Thome! d 1987b . with the remnant mass of ~ 100M Q . 

For comparison, in Figure [2] we also plot the spectrum as- 
sumed in previous w orks (dot-dashed; Buo nanno et alj[2005b 
ISandick et al. 2006). The result of our calculation is very dif- 
ferent from their spectrum. This is due to consideration of hy- 
drodynamical features. Pop III stars have higher temperature 
prior to collapse, which leads core-bounce by the (gradient of 
the) thermal pressur e and not by the nuclear forces as in or- 
dinary supernovae dFrver et alJl2001l:ISuwa et alj|2007l) . Such 
thermal bounce makes the central density lower at the time of 
bounce (<9(10 12 )g cm" 3 ) so that the dynamical timescale (oc 
p~ 1//2 ) becomes longer, leading to the smaller time-variability 
in the waveforms and the smaller typical frequencies (a few 
10 Hz) where the matter contribution peaks, which are in the 
1 kHz regime for ordinary supernovae. 

It should be noted that current numerical simulations only 
encompass a few seconds at most and do not cover the strain 
spectrum below a fraction of a Hertz. For such a l ow fre- 
quency region, we apply the zero-frequency limit dSmarrJ 
119771: lEpsteinl 1 1978J) to e xtract the GW waveforms as done 
in [Buonanno et al. ( 2005b- 

Finally, we discuss the GW emission from stars with differ- 
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FIG. 2. — Spectral distributions of the gravitational waves from anisotropic 
neutrino emissions (blue dotted line), matter motions (green dashed line), 
and total (red solid line). Due to phase cancellations between the matter and 
neutrino components, the amplitude of the combined spectrum is smaller than 
the dominant components in certain frequency regions. The dominance of the 
neutrino-oriented gravitational wave component at frequencies below ~ 10 
Hz is clearly seen. The dot-dashed line represents the spe ctrum employed in 
previous works iBuonanno et al. 2005; Sandick et al. 2006). 

ent mass models. It is noted that the qualitative features are 
common to all models. In Table[T] we show the characteristic 
values of GW spectra. It can be seen that the peak amplitude, 
hmax, increases with the mass of the star, while the peak fre- 
quency, / m ax, remains roughly the same. The frequency where 
the neutrino and matter contributions become equal,/ eq , de- 
creases with increasing mass. The GW amplitudes for the 
zero frequency limit, Df\h\(f — > 0), increases with the mass 
since the luminosity of neutrinos becomes larger. 



TABLE 1 

Characteristic Quantities of Gravitational Wave 



Mass 




/max 


/eq 


Df\h\(f->0) 


(M ) 


(km) 


(Hz) 


(Hz) 


(km) 


300 


0.66 


21 


7.8 


0.30 


500 


0.88 


19 


6.0 


0.46 


700 


0.86 


23 


4.4 


0.38 


900 


0.99 


18 


4.0 


0.54 



4. GRAVITATIONAL WAVE BACKGROUND FROM POP III STARS 

We are now in a position to discuss the contribution of GWs 
from Pop III st ars to the background radiation. According to 
lPhinnevl ( l20 01). the sum of the energy densities radiated by a 
large number of independent Pop III stars at each redshift is 
given by the density parameter flcw(f) = P~ c l (d PGw/d log f) 



m / dm0(m)Z) 2 / ,3 |H/')| 2 



l5Gp c J l+z 

(2) 

where p c is the critical density (3// ( 2 /(87rG)), ip(z) is the 
star formation rate (SFR), (j)(m) is the initial mass function 
(IMF) of Pop III stars, and /' is t he red shifted freque ncy, 
(l+z)/- We employ model 2b of ISandick etail (120061) for 
the SFR, which is appropriate for very massive stars, from 
270M Q to 500M q. As for the IMF, w e employ the same pa- 
rameterization as Sandic k~et al.l (120061) . 4>(m) oc m~ 23 , which 
is normalized by J dm<fi(m) = 1. The cosmological model 
enters with \dt/dz\ = [(1 + z)H(z)]~ l and, for a flat geome- 
try, H(z) = H [Q\ + Sl m (l + z) 3 ] 1/2 . We use the parameters 




FIG. 3. — The energy density parameter of gravitational wave background 
(GWB). The upper line of the shaded region is the result of our calculation. 
The lower line is for the renormalized SFR of//// ~ 10~ 5 . In low frequency 
region (/ ^ 0.1Hz), f!ow / because f\h\ is constant (see Figure 2) and 
the integrand of Eq. (2) is / 3 \h\ 2 oc / in this region. The horizontal dashed 
line shows the GWB produced during slow-roll inflation assuming T/S = 0.3 
for the ratio of the tensorial and scalar contributions to the cosmic microwave 
background radiation anisotropy and no running of t he tens orial power-law 
index, which is evaluated from Eq. (6) of ITurnei (19971) . The sensitiv- 
ity curves of space - an d (correlated) grou nd-based detectors are taken from 
IKudoh et all <2006T) and lBuonannTj <2003l) . 

ft A = 0.7 3 A = 0.27, and H p = 100 h km s" 1 Mpc" 1 with 
h = 0.7 1 (Spe reelet al1l2007h . 

In Figure [3] the calculated f2cw is plotted with the sensi- 
tivity curves of future detectors. The upper edge of the red 
shaded region is obta ined by our calculation with model 2b of 
Sandi cket alj (|2006). This corresponds to almost the upper 
limit of the SFR because the stars with the mass range con- 
sidered in this model entirely collapse to BHs and do not con- 
tribute to the chemical evolution of their environment. In this 
way, a high SFR is obtained without metal overproduction. 
Meanwhile, the lower limit is the same shifted downwards by 
a factor 7000, which corresponds to the baryon fraction of Pop 
III stars of 10~ 5 , suggested by Pop III star formation theory 
(Omukai, private communication). It can be seen that the con- 
tributions from neutrinos, which dominates in low frequency 
region (~ 1 Hz), are within the detection limit of the planned 
detectors DECIGO and BBO. In contrast, detection by LIGO 
III, which is a ground based detector, may prove difficult be- 
cause the amplitude in the high frequency region is smaller 
than the current estimated detection limit. Even if the SFR of 
Pop III stars is as small as in the case of a baryon fraction of 
~ 10~ 5 , the GWB from Pop III stars might be within the de- 
tection limit of ultimate-DECIGO. The GWB spectrum pre- 
dicted in this letter is larger than previous works b y a factor of 
40 in the low frequency region (see Figure 8 of Sandick et al. 
20061 which represents larger amplitude than lBuonanno et al] 
2005). This is due to the difference of a GW spectrum of a 
single Pop III star collapse as already mentioned in §3. De- 
pending on the baryon fraction of Pop III stars, the GWB 
from Pop III stars might give a strong contribution, masking 
the GWB generated in the inflationary epoch (the horizontal 
dashed line). We furthermore point out that the amplitude of 
the GWB is highly dependent on the SFR but less sensitive to 
the index of IMF because the GW emission of different mass 
Pop III stars is similar in the mass range focused on in this 
letter. 

5. SUMMARY AND DISCUSSION 

In this letter, we have presented the GWB spectrum from 
Pop III stars based on the hydrodynamic core-collapse sim- 
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ulations of Pop III stars (ISuwa et al J 12007). Extracting the 
GWs contributed from mass motions and anisotropic neutrino 
radiation, we have found that the GW emission from neutri- 
nos is dominant over the ones from matter at low frequencies. 
Using the spectrum of single Pop III stars, we have calculated 
the GWB spectrum by summing up the contribution from in- 
dividual stars. We have shown that the amplitudes from Pop 
III stars might be large enough to be detected by interferome- 
ters such as DECIGO and BBO. 

The predicted GWB is in the sensitivity range of future 
planned detectors so that the Pop III stars might limit the sen- 
sitivity of such detectors. This is analogous to the discus- 
sion about the the Galactic population of close white dwarf 
binaries, which could produce a confusion-limited GW signal 
at the lower end of the sensitivity band of LISA dHils et alj 
fl990l) . 

It is noted th at the Newtonian simulation in the paper 
dSuwa et al.l 120071) is nothing but an idealized study to de- 
scribe the dynamics with BH formation. It is by no means 
definitive. To explore these phenomena in more detail will re- 
quire fully general relativistic simulations and much more bet- 
ter neutrino transfer. However, in comparison with the forgo- 
ing GWB studies simply using the GW waveform of ordinary 
core-collapse supernovae as a template, we have found sig- 
nificant differences in the properties of the waveforms, which 
made it possible to realize the potential importance of Pop III 
stars as a contribution to the GWB. As mentioned, our results 
show enhancements in the GWB amplitudes from 0.1 to 1 Hz 
compared to earlier works, where the currently planned detec- 
tors are most sensitive. 

Finally, we shall briefly discuss the uncertainties of our 
models. Since little is known about the angular momentum 
distributions of Pop III stars, we took the initial rot ation rate of 
extrem ely rapidly rotating cores of massive stars lHeger et all 
(2000) as a reference suggested by the study that Pop III stars 
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could rotate rapidly due to the insuffi cient mass-loss in the 
main sequence stage (He ger et alJ20 03). Obviously more sys- 
tematic studies changing the initial rotation rates and profiles 
are required. Next, the IMF of Pop III stars is also uncertain. 
It is true that the mass range employed here (300- 1000M©) is 
preferred by recent studies of Pop III star formation, but that 
the possible modification of the IMF leads t o large changes 
of the predicated GWB s pectrum (see, e.g, ISchneider et alj 
l2000Hde Araujo et alJ2002h . This means that the GWB could 
be a powerful tool to investigate the IMF. Although the baryon 
fraction of Pop III stars taken here is indeed not too large for 
the explanation of the infrare d background exce s s by the UV 
photo ns from Pop III stars ( Sant os et alJ 120021 : iDwek et all 
120051) . there still remains large uncertainty. The last uncer- 
tainty is the redshift where the SFR becomes maximum. It 
should be mentioned that the frequencies where the GWB 
from Pop III stars peaks are sharply dependent on this red- 
shift. Better correlation analysis between multiple spaced- 
based interferometers should be very helpful, because this 
can im prove the sensitivity of GWB detection (Kudoh et al. 
2006). All these results suggest that detections of GW back- 
ground from Pop III stars can be an important tool to supply 
information about the formation history of Pop III stars. 
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